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As a generalization of the Baxter-Wu model, we numerically study the Ising model with triplet 
interactions on the Union- Jack lattice, which we call as Union-Jack Baxter-Wu model. By means 
of finite-size scaling analysis based on Monte Carlo simulations and transfer matrix calculations, we 
numerically determine the critical exponents (or scaling dimensions) and the central charge of the 
model. The critical exponents studied include y t , yn, Vhi, and where yt governs the critical 
behavior of the correlation length, yu governs the critical behavior of the magnetization of the 
whole lattice, yni and yui govern respectively the critical behaviors of the magnetizations on the 
two sublattices of the Union- Jack lattice. For the critical exponents y t , yni, Vh2, and the critical 
points of the model, our numerical estimations coincide with the exact solutions. For the critical 
exponent y^ and the central charge c, as we know, there is no exact solutions, thus our numerical 
results are the first determination of these exponents. The exact solutions and our numerical results 
show that the critical exponent yhi and the central charge c take the same values as those of the 
4-state Potts model, which suggests that the model 'partially' conserves the critical properties of 
the 4-state Potts universality. 

PACS numbers: 05.50.+q, 64.60.Cn, 64.60.Fr, 75.10.Hk 
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I. INTRODUCTION 



The exact solution|l|, |2[ of the Ising model significantly promotes the research of phase transition and critical 
phenomena. After that, the Ising model becomes one of the most famous lattice models in statistical physics. 
Another famous Ising system, known as the Baxter- Wu modelQ, is defined on the triangular lattice. The model 
has three-spin interactions, and the reduced Hamiltonian is written as 

-H/k B T = K Wfc> (!) 

<i,j,k> 

where Si, Sj, Sk are the Ising spins sitting at the three points of an elementary triangle, and K is the coupling 
constant of the three spins. The sum takes over all the elementary triangles. This model was firstly proposed by 
Wood and GriffthsQl and exactly solved by Baxter and Wu0] by relating the model with the coloring problem on 
the honeycomb lattice. The solution gives the critical exponents yt = 3/2 (a = 2/3) and yh = 15/8 (i] = 1/4), which 
are exactly the same as those of the 4-state Potts model[l, @, which can be derived by the Coulomb gas theoryQ. 
This means the Baxter- Wu model is in the same universality of the 4-state Potts model. This result can also be 
understood by the fact that both the Baxter- Wu model and the 4-state Potts model have four-fold degenerate ground 
state. The critical properties of the 4-state Potts model are affected by the logarithmic corrections [§| because of 
the second temperature field, which are marginally relevant, however, this is not the case for the Baxter- Wu model. 
Deng et. al. generalize the Baxter- Wu model in [9j where the spins are allowed to be q states (q can be larger than 
2) and the up- and down-triangles can have different coupling constants. Both generalizations lead to discontinuous 
phase transitions. 




FIG. 1. The Baxter- Wu model on the Union- Jack lattice: the lattice consists of two sublattices, the spin on sublattice A is 
denoted as a and the spin on sublattice B is denoted as s; for a square unit cell, there are four spins 0-1,0-2,03, and 0-4 on 
sublattice A and one spin s on sublattice B. 



In 1972, Hintermann and Merlinifioj considered an Ising system with pure three-spin interactions on the Union- 
Jack lattice (as shown in Fig. [1]) 

- H/k B T = 2J s(Kiaia 2 + K^cj^z + K^a^ + K^a^ai), (2) 

{s,a} 

where the sum takes over all the square unit cells and Kt (i — 1, 2, 3, 4) are the coupling constants. The model can 
also be viewed as a generalization of the (triangle) Baxter- Wu model, thus in current paper, we call the model as 
Union- Jack Baxter- Wu model. For the ferromagnetic case [Ki > 0), the model also has four-fold degeneracy of the 
ground state (as shown in Fig. [5]), thus 4-state Potts universality is expected for the model. However, the exact 
solutions [lol show that the model has critical exponents which change along with the coupling constants. 

Through mapping to the eight- vertex model [HI, [u|, this model can be exactly solved 10], the critical point and 
critical exponent y t are obtained. Using the similar method, Wu calculates the spontaneous magnetizations of the 



model on the two sublattices, the results show that the two magnetizations possess different critical exponentsjll|. 

In current paper, we will numerically study the ferromagnetic Union- Jack Baxter- Wu model for the case with 
K\ = K3 > and K% = K4 > 0, the numerical procedure includes the Monte Carlo simulations and the transfer 
matrix calculations. The paper is arranged as follows: In Sec. [TTl we summarize the results of exact solutions. In 
Sec. IIIH we describe the algorithm that we used in Monte Carlo simulations and give the associated numerical 
results. In Sec. IIV1 we introduce the transfer matrix calculations and give the associated numerical results. In Sec. 
fVl we summarize the paper and give some discussions. 
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FIG. 2. (Color online) Four ground states of the Union-Jack Baxter- Wu model. 



II. EXACT SOLUTIONS 

We summarize the results of exact solutions in this section [T3, EH- 

Assign arrows between the nearest-neighboring Ising spins of sublattice A (Fig. [3]) : if the two spins besides the 
arrow is the same, the arrow is rightward (or upward), otherwise it is leftward (or downward). There is a two-to-one 
correspondence between the Ising spin configurations and the arrow configurations. 
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FIG. 3. (Color online) The mapping of the Union-Jack Baxter- Wu model to the eight-vertex model and the Boltzmann 
weights of the vertexes. The Ising spins shown here are the spins on sublattice A, namely the a spins. 



After taking sum over the spin s (in the center of the unit cell) , the Boltzmann weights of the vertexes are 

a = cosh(A"i +K 2 +K 3 + K 4 ), 
b = cosh^i - K 2 + K 3 - Ki), 
c = cosh(A^i -K 2 -K 3 + K A ), 
d = cosh(A^i + K 2 - K 3 - Ki). 



(3) 



Thus the model is mapped to the symmetric eight vertex mo del [ 13. Il3j . which was exactly solved by Baxter. For 
the ferromagnetic case that we consider in current paper, the critical point is given by 



a = b + c + d, 

which gives K c — log(l + V2)/2 for the uniform case K = K\ = K 2 = K 3 = Ki. 



(4) 
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The singularity of the free energy density is governed by 

f sing oc\T-T c \^ u = \T-T c \ 2 /y* (5) 

when 7t/2m is not integer (The case with ir/2u = integer is not considered in current paper). T c = 1/K C is the 
critical temperature and < u < ir is given by [Tl| 

cos u — — tanh[- log — ] if a > b, c, d, or b > a, c, d (6) 

= tanh[- log — ] if d > a, b, c, or c > a, b, d. (7) 

For the ferromagnetic case that we consider, u is determined by ©• © gives the critical exponent 

y t = 2 A (8) 

where yt also governs the critical behavior of the correlation length: £ ~ \T — T c \~ x l yt . For the uniform case, (JSJ) 
gives y t = 4/3. 

A very interesting critical property of the Union- Jack Baxter- Wu model is: the spontaneous magnetizations of 
the two sublattices possess different critical exponents 

M A <x{T c -Tf«, (9) 
M s oc (T c - Tf >, (10) 

where Ma and Mb are the magnetizations on sublattice A and B, respectively. The critical exponents f3 a and f3b 
are obtained by Wu[li1] 

Pa ~ ~rp. j 
16u 

Notation: in Eq. (16) of [ll|, the author reversely writes the critical exponents of Ma and Mb- 
According to the scaling law, the magnetic exponents yu\ and yu2 are 

15 

Vhl = 2 - PaVt= "g-j 

y, 2 =2-/3 b y t =^±^. (12) 

Here y^i has fixed value 15/8, while yhi varies continuously with the parameters of the model. In the following 
sections, we will numerically study the critical properties of the Union- Jack Baxter- Wu model, the critical points 
and critical exponents will be numerically verified. Especially, Wu's results of /? and /3b (consequently yhi and y^i) 
are based on the conjectured results of the spontaneous magnetizationflij and spontaneous polarization [15] of the 
eight- vertex model, thus numerical verification is very necessary. 

We are also interested in the critical property of the magnetization of the whole lattice M, whose critical behavior 
is assume to be governed by a critical exponent j3 (consequently yh = 2 — (3y t ) 

M cx (T c - Tf. (13) 

As we know, the critical exponent yh (or /?) is still neither exactly solved nor numerically determined, we will 
determine it by Monte Carlo simulations in the next section. 



III. MONTE CARLO SIMULATIONS 

The triangle Baxter- Wu model has been simulated by the Metropolis algorithm [l6j. the Wang Landau algorithm [l7j. 
and a cluster algorithm (9l [l8| . The cluster algorithm Q is similar to the Swendsen-Wang algorithm fl9l| for the Potts 
model. In this algorithm, the triangular lattice is divided as three triangular sublattices. By randomly freezing the 
spins on one of the sublattices, the other two sublattices compose a honeycomb lattice with pair interactions, then 
a Swendsen-Wang type algorithm can be formed to update the spins on this honeycomb lattice. For details of this 
cluster algorithm, see js}. In current paper, we suitably modify this cluster algorithm to simulate the Union- Jack 
Baxter- Wu model. 

The algorithm is divided into three steps: 
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1. Step 1, update the spins on sublattice A. 

Let the spins on sublattice B be unchanged (frozen) in this step, the interactions of the spins on sublattice A 
reduce to pair interactions. 

A. bonds. Every edge on sublattice A is clamped by two triangles, the interactions of the left and right 
triangles can be denoted respectively as Ki and K r , where Ki and K r may take the values Ki, K2, K3, and 
K4. Place a bond on this edge with probability p = 1 — e~ 2Kl if only the product of the three spins of the 
left triangle is 1, p = 1 — er 2Kr if only the product of the right triangle is 1, p = 1 — e -2K t -2K t ^ both the 
product of the left triangle and the product of the right triangle are 1, p = otherwise. 

B. clusters. A cluster is defined as a group of sites on sublattice A connected through the bonds. 

C. update spins. Independently flip the spins of each cluster with probability 1/2. 

2. Step 2, update the spins on sublattice B and half of the spins on sublattice A. 

This step is very similar to step 1, but we freeze only half of the spins on sublattice A, which are labeled 
Ai (Fig. 2]). The other spins (on sublattice B+A 2 ) are updated, whose interactions also reduce to pair 
interactions when Ai spins are frozen. 

3. Step 3, update the spins on sublattice B and the other half of the spins on sublattice A. 

In this step, the spins labeled A2 (in Fig. [4]) are frozen , other spins (on sublattice B+Ai) are updated. 

In a complete sweep, all the spins are updated twice. 

A 2 A, A 2 Ai 

"XlXlXJ 
|X|X|XI 

FIG. 4. (Color online) When B spins are frozen, the interactions of A spins (Ai and A2 spins) reduce to pair interactions; 
when Ai spins are frozen, the interactions of B + A2 spins reduce to pair interactions; when A2 spins are frozen, the 
interactions of B + Ai spins reduce to pair interactions. 

In the simulations, the sampled variables include the specific heat C, the magnetization of the whole lattice M, 
the magnetization on sublattice A Ma, and the magnetization on sublattice B Mb- The specific heat is defined by 
the fluctuation of internal energy 

C= J_ <*>-<*>' (14) 
L 2 k B T 2 ' { ' 

where E is the internal energy, L is the linear size of the system, and (■ • • ) means the ensemble average. The 
magnetizations are defined as 

N N 

(I E °i + E s *\) 

M= 1=1 m * =1 , (15) 

iV 

<|£*.|> 

Ma = , (16) 

JV 

<|E^|) 

Mb = ^ , (17) 



G 



where N = L 2 is the total sites on sublattice A or B. 

In order to demonstrate our numerical procedure based on the Monte Carlo simulations, we take the uniform 
case (K = K\ = K 2 = K3 = K4) as an example. The cluster algorithm is very efficient, it easily allows us 
to do meaningful simulations for system with linear size up to L — 256. All the simulations are performed at 
the critical point K c . After equilibrating the system, 10 7 samples were taken for each system size. Using the 
Levenberg-Marquardt algorithm, we fit the data of C according to the finite-size scaling formula (20T| 



C(L) = C + L 2yt - d (a + bL yi ), 



(18) 



where bL Vl are the leading correction-to-scaling and y\ < is the leading irrelevant exponent, d = 2 is the dimension 
of the lattice. Co, a, and b are unknown parameters. The fitting yields y t — 1.332(2), which is in good agreement 
with the exact solution yt = 4/3. 

Similar procedure was applied to the magnetizations M, Ma, and Mg, according to the finite-size scaling 
formulas 12(1 



M = L Vh - d {a-\ 
Ma = L yhl - d {a- 
M B = L Vh2 - d (a- 



-bL Vl ). 



(19) 
(20) 
(21) 



The log-log plot of M, Ma, and Mb versus L are shown in Fig. [5] The fitting of the data yield yh — 1.857(1), 
yhi — 1.874(1), and yh2 = 1.833(1), where yhi and yya are in good agreement with the exact solutions yhi = 15/8 
and yhi = H/6, respectively. 



0.55 




192 256 



FIG. 5. (Color online) Log-Log plot of the magnetizations (M, Mi, M2) versus system size L for the Union- Jack Baxter- Wu 
model at the critical point K c = log(l + v / 2)/2 for the uniform case (K — K\ = K2 = K3 = K4). 

We also simulate the cases with r = K' /K =1, 2, 3, 4, and 5, where K = K\ = K3 and K' = K2 = K4. All 
the exact solutions and numerical estimations are listed in Table HI From the table we can see that the numerical 
estimations for y t , yhi, and yhi are in good agreement with the exact solutions. To the best of our knowledge, there 
is no exact solution for yh until now, thus our numerical result is the first determination of this critical exponent 
for the Union- Jack Baxter- Wu model. 



IV. TRANSFER MATRIX CALCULATIONS 



As shown in Fig. [51 we define the row-to-row transfer matrix 



Tg,gi = ]^[exp[si(ificri<Ti+i + K 2 a l+ iCT l i+l + K3a' i+1 a' i + K i n' i a l )] 

{si,« 2 ,— «i} i=1 
L 

= 2 L ^Qcosh(i4Ticr i cr i+ i + K 2 (Ti+i(j' l+1 + ir 3 cr- +1 (T^ + K 4 a' l a i ), 
i=i 



(22) 
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TABLE I. Critical points, critical exponents, and scaling dimensions of the Union- Jack Baxter- Wu model, r = K' /K, 
where K — K\ = Kz and K 1 = if 2 = Ka- E=Exact result, MC=Numerical result based on Monte Carlo simulations. 
TM=Numerical result based on transfer matrix calculations. The critical exponents and the scaling dimensions are related 
by y t = 2 - X t , y h \ = 2 - X hl , and y h2 = 2 - X h2 . 



V 


1 


2 


3 


4 


5 


K c E 
TM 


0.4406867935 
0.44068679(5) 


0.3046889317 
0.3046889(1) 


0.2406059125 
0.2406059(1) 


0.2017629641 
0.2017629(1) 


0.1751991102 
0.1751991(2) 


yt e 

MC 

X t E 
TM 


4/3 

1 332(2) 
2/3 

0.66666(1) 


1.39668184 
1.397(4) 
0.60331816 
0.603318(1) 


1.47604048 
1.478(4) 
0.52395952 
0.52396(1) 


1.53960311 
1 540(31 
0.46039689 
0.46040(2) 


1.58921160 
1 590(4) 
0.41078840 
0.41079(1) 


Vhi E 
MC 

X hl E 
TM 


15/8 

1.874(1) 

1/8 

0.125000(1) 


15/8 

1.874(2) 

1/8 

0.125000(1) 


15/8 

1.876(2) 

1/8 

0.125000(1) 


15/8 

1.877(3) 

1/8 

0.12500(1) 


15/8 

1.873(4) 

1/8 

0.12498(3) 


2fo2 E 

MC 

X h2 E 
TM 


11/6 

1.833(1) 

1/6 

0.166666(1) 


1.84917046 
1.848(2) 
0.15082954 
0.150829(1) 


1.86901012 
1.870(2) 
0.13098988 
0.130990(1) 


1.88490078 
1.886(3) 
0.11509922 
0.115099(1) 


1.89730290 
1.896(4) 
0.10269710 
0.102696(2) 


W. E 
MC 


1.855(2) 


1.862(3) 


1.874(3) 


1.882(4) 


1.888(4) 


c E 
MC 


1 

1.000000(1) 


1 

1.000000(1) 


1 

0.99998(3) 


1 

1.0000(2) 


1 

0.9999(1) 




FIG. 6. (Color online) Definition of the row-to-row transfer matrix, where periodic boundary conditions are applied, namely 



where a = (ci , o"2 ) • ' " ,<Jl) and 3' — (cr^er^''' are the states of two neighboring rows. In the definition, 

periodic boundary conditions are applied, namely <Tl+i — <J\ and cr' L+1 — a[. For a system with M rows, the 
partition sum is 



7 1 - - T- - ■ ■ ■ T- 



T~ - T- - • • • T- - 

f 1 ,f 2 ""2 3 fM)fM+l 



^ ] ■'■ ' T i .' T :; '"j;."! 



= Trace(T M ), 



(23) 



where CTj is the state of the i-th row and periodic boundary conditions are also applied, namely (Jm+i = In the 
limit M — > 00, the free energy density is given by 

1 



/= nm -^-logTrace(T M ) 



1 



= ton — log(Af +A£ 



ML 
logAi 



A 



(24) 



where Ai, A2, A3, • • • are the eigenvalues of T, with Ai > A2 > A3 • • • . An exact solution of Ai determines the free 
energy density of the system. Furthermore, not only the free energy density, but also some critical exponents are 
also determined by the eigenvalues of the transfer matrix, which will be shown concretely bellow. 
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To numerically calculate some of the eigenvalues of the matrix, the sparse matrix technique is very useful, it 
sharply reduce the requirement of computer memory for storing the matrix elements. For detailed study of this 
technique, see (2l| - [23 |. For the Union- Jack Baxter- Wu model, the dimension of the matrix T is — 2 L . By using 
the technique of sparse matrix decomposition, we are able to calculate the eigenvalues of T with L up to 22. In 
order to obtain meaningful eigenvalues, the L must be even because of the symmetry of the ordered phase, as shown 
by Fig. [2] especially by Fig. [2jc) and (d). For a finite system, the symmetry can be conserved only if the size is 
even. 

In order to study the critical properties of the Union- Jack Baxter- Wu model, we calculate the scaled gap X(K, L), 
which is generally defined as 

X(K,L) = ±-log(^), (25) 

where Ai and Aj are the leading and the other eigenvalues of the matrix T, respectively. According to the conformal 
invariance theory[25j, the scaled gap X(K,L) is related with the correlation length £(K,L) with 

where £(K, L) governs the decay of the correlation function G(r): 

G(r) oc e - r '^ K ' L \ (27) 

The correlation function G(r) can be energy correlation or magnetic correlation. For example, the magnetic cor- 
relation function can be defined as G(r) = (sos r ) or the probability of two spins with distance r have the same 
value. For the Union- Jack Baxter- Wu model, because we arc interested in the different critical properties of the 
two sublattices, we can define other types of magnetic correlation functions, where we confine sq and s r on one of 
the sublattices. For clarity, we let G(r) represent the magnetic correlation function on the whole lattice, let Ga{t) 
and Gb{t) represent respectively the magnetic correlation functions on the A and B sublatices. We further express 
the scaled gaps corresponding to G(r), GA{r), Gs(r) as Xh(K, L), Xhi(K, L), Xhi{K, L), respectively, and express 
the scaled gap corresponding to the energy correlation function as X t (K, L). 

In the procedure of transfer matrix calculations, we find that, the scaled gaps are determined by the first four 
eigenvalues as (f25|) . The concrete order are determined by the magnitude of the gaps. For example, when Xhi < 
Xhi < Xt, the scaled gap Xhi(K, L) is determined by eigenvalues Ai and A 2 , Xh%{K, L) is determined by Ai and 
A3, and X t (K, L) is determined by Ai and A 4 : 



L 
2^ 


log 


A x 
A 2 


L 


log 


Ar 


2^ 


A 3 


L 
2^ 


log 


Ar 

A 4 


by 






L 


log 


A x 


2^ 


A 3 


L 
2^ 


log 


A x 
A 2 


L 


log 


Ai 


2^ 


A 4 



(28) 



X hl (K,L) 



X t (K,L) = —lo gl ^, (29) 
and so forth. 

Unfortunately, we are unable to find which eigenvalues determine Xh(K,L). 

According to the finite-size scaling[2(| and the conformal invariance [25j| . the scaled gaps Xhi and Xh2 in the 
vicinity of the critical point K c satisfy 

X hl (K, L) =X hl +a(K- K c )L yt + buL^ + ■■■, (30) 
X h2 {K, L) = X h2 + a'{K - K e )V* + b'uL^ + ■■■ , (31) 

where Xhi and Xh2 are the scaling dimensions, u is the irrelevant field, and y u < is the associated irrelevant 
exponent, a, a', b, and b' are unknown constants. Fig. [7] and Fig. [S]are illustrations of Xh^{K, L) and Xh2(K,L) 
versus K respectively, for the uniform case. 
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FIG. 7. (Color online) Scaled gap Xhi(K,L) versus K for a sequence of system size L for the uniform Union-Jack Baxter 
Wu model (r = 1), whose critical point is K c = log(l + y/2)/2 = 0.440687. 
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FIG. 8. (Color online) Scaled gap Xh2(K,L) versus K for a sequence of system size L for the uniform Union- Jack Baxter 
Wu model (r = 1), whose critical point is K c = log(l + \/2)/2 = 0.440687. 



We numerically solve the finite-size scaling equation 

X hl {K,L)=X hl [K,L-l), 



(32) 



X h2 (K,L)=X h2 (K,L-l), 
and denote the solution by K C (L), which satisfies 

K C (L) =K C + auL v "- yt + ■■■ 



(33) 



(34) 



Here a are unknown constant. y u < and y t > 0, thus K C (L) converges to the critical point K c for increasing 
system size. (|34p is used to determine the critical point in our numerical procedure. For the uniform case (r = 1), 
it gives K c = 0.44068679(5), which is in good agreement with the exact solution K c — log(l + v / 2)/2. We have also 
estimated the critical points of the cases with r = 2,3,4, and 5, all of our numerical estimations of K c , which are 
consistent with the exact solutions in a high accuracy, are listed in Table HI Sincerely, both (|30)) and (|3T|) can be 
used to estimate the critical points, in Table U we list the best estimations. 
Exactly at the critical point K c , (f30|) and (|3l"j) reduces to 

(35) 
(36) 



X hl (L) = x hl + buLy- 

X h2 (L) = X h2 + b'uLV« 



These equations are used to determine the scaling dimensions X^i and Xh 2 in our numerical procedure. For the 
uniform case, it gives Xhi = 0.125000(1), and Xh 2 = 0.166666(1), which are consistent with the exact solutions 
(Table |T|. We also applied the same procedure to the other cases, the numerical estimations and exact solutions are 
listed in Table Q] 
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FIG. 9. (Color online) Scaled gap X t (K, L) versus K for a sequence of system size L for the uniform Union- Jack Baxter Wu 
model (r = 1), whose critical point is K c = log(l + y/2)/2 = 0.440687. 



We also calculate the scaled gap X t (K,L), its critical behavior seems somewhat different from that of Xhi or 
Xh2- This is shown in Fig. [H] for the uniform case. The lowest point of the line for a given L can be denoted as 
(K C (L), Xt(K c (L), L)), and the value of K C (L) seems to converge to the critical point K c for increasing L. This 
property can also be used to estimate K c , for examples, see [2o| . But we haven't do this in current paper, because 
we have estimated K c by the critical behavior of Xhi or Xya- However, exactly at the critical point K c , the value 
of X t (K c ,L) seems to satisfy the formula 

X t (K c ,L) = X t + buLV~ + --> , (37) 

which is very similar to (|3"5)) or (|3T>1) . We use (|3T|) to estimate the scaling dimension X t in our numerical procedure, 
the results are also listed in Table [I] 

The scaling dimensions Xhi,Xh2, and X t are related with the critical exponents bv[2o| 





= 2 






= 2 


— Xh2 


Vi 


= 2 


-x t . 



(38) 

This allows us to compare the results of the transfer matrix calculations, the Monte Carlo simulations, and the 
exact solutions (Table fl}. 

According to the finite-size scaling (2^1 and the conformal invariance(25|. the free energy density scales as 

f(L) ~ /(oo) + ^ (39) 

at the critical point, where c is the central charge in the conformal theory (25|. Fitting the data of the free energy 
density according to (I3T))) , we get the central charge c = 1.000000(1) for the uniform case. For the other cases, the 
numerical estimations of the central charge are also listed in Tabic |TJ As we know, this is the first determination 
of the central charge of the Union- Jack Baxter- Wu model. It seems that the central charge takes the fixed value 
c = 1, which is the same as that of the 4-state Potts model. 



V. CONCLUSION AND DISCUSSIONS 



In conclusion, we have numerically studied the critical properties of the Union- Jack Baxter- Wu model, by means 
of finite-size scaling analysis based on the Monte Carlo simulations and the transfer matrix calculations. For the 
critical points and critical exponents m ( Xt), yui (Xhi), and y%2 (Xya), our numerical estimations are in good 
agreement with the exact solutions [Tol I 111]. To the best of our knowledge, there is no exact solutions of the critical 
exponents yh and the central charge c, thus our numerical results are the first determination of these exponents. 
The numerical results show that yh varies continuously with the parameters of the model, while the central charge 
seems to take the fixed value c = 1, which is the same as that of the 4-state Potts model. 

The ground state of the (ferromagnetic) Union- Jack Baxter- Wu model is four-fold degenerate, thus the phase 
transition of the model is expected to be in the universality of the 4-state Potts model. However, the model has 
much more complicated critical properties than the 4-state Potts model, which are shown by the continuous change 
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of the critical exponents {yt, yh, and y/12), and by the difference of the critical exponents on the two sublattices (yhi 
and y h2 ). 

It seems that the model 'partially' conserves the critical properties of 4-state Potts universality, which are shown 
by the critical exponents y^i and the central charge c, whose values are fixed and the same as those of the 4-state 
Potts model. 

The critical properties of the 4-state Potts model are affected by the logarithmic corrections [8j, however, in our 
numerical procedure, all the data are well fit without the logarithmic correction. 
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